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Abstract 

We study the role of the A-isobar current on the response functions for high energy inclusive 
quasielastic electron scattering from nuclei. We consider a general Lagrangian which is compatible 
with contact invariance and perform a fully relativistic calculation in first-order perturbation the- 
ory for one-particle emission. The dependence of the responses upon off-shell parametrizations is 
analyzed and found to be mild. A discussion of scaling behaviour and a comparison with various 
non-relativistic approaches are also presented. 
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1 Introduction 



In recent work [1, 2, 3] we have studied the role of pionic correlations for electron- nucleus 
scattering in the kinematical domain of the quasielastic peak, focussing our attention on the 
lp-lh sector. The calculation has been carried out on the basis of a Relativistic Fermi Gas 
(RFG) model in which an exact relativistic treatment of the problem can be accomplished. 
In the domain of perturbation theory, all diagrams including one pionic line were included 

- in particular the meson-exchange currents (MEC) of contact and pion-in-flight types - 
so that the corresponding current was automatically gauge invariant. 

In [1, 2] the diagrams involving the excitation of a virtual A-resonance, represented in 
fig. 1, were neglected. Although they do not affect the gauge invariance of the theory, since 
they are associated with a conserved current, it is well-known that such contributions do 
modify the nuclear response functions, especially in the transverse channel. The aim of this 
paper is to extend the model of refs. [1, 2] to include the A-induced meson-exchange-currents 
of fig. 1, in order to have a complete understanding of pionic effects in the quasielastic 
peak (QEP) as far as the lp-lh channel is concerned. In contrast with non-relativistic 
calculations [4] our results do not involve any expansions in energy-moment um/M, where 
M is a typical baryonic mass, and can therefore be applied when studying the response at 
high momentum and energy transfers. 

The present relativistic treatment of the A current allows one to study several aspects of 
A electroexcitation which are of special interest theoretically. In particular, in the effective 
Lagrangian approach to electroexcitation of the A resonance it is known that there is free- 
dom at the electromagnetic (EM) vertex due to the off-shell behavior of this resonance [5]. 
In contrast to earlier approaches to the A-exchange current in a relativistic model [6], here 
we consider a more general '-/NA interaction Lagrangian. First of all, current conservation 
restricts the form of the vertex to a superposition of three covariants, for which one can 
use any of the choices described for instance by Jones and Scadron [7]. One possible choice 
is the familiar set of magnetic dipole, electric quadrupole, and Coulomb quadrupole multi- 
poles, used in the pion electroproduction analyses of [8, 9]. In this work we use instead the 
standard "normal parity" set, analogous to the Dirac-Pauli decomposition of the nucleon 
form factor [10]. The corresponding term in the traditional chiral Lagrangian of Peccei [11] 
is just a particular case of the normal parity set with two of the three terms equal to zero. 

Going a step beyond the Jones-Scadron vertex, the off-shell propagation of massive vector 
fields with spin S > \ has been exhaustively discussed in the literature [5, 12, 13]. In 
particular, a spin | field generates contributions involving the S —\ sector of this field in the 
effective amplitudes. In the case of the A, this is obtained by constructing the most general 
vertex which is invariant under a special contact symmetry of the free Lagrangian (see [5] for 
details). The invariance of the 7./VA and 7riVA vertices under this "point transformation" of 
the field requires the introduction of additional parameters in the Lagrangian. Attempts have 
been made to fix some of these "off-shell" parameters by fitting the pion electro- and photo- 
production data [14] and, more recently, Compton scattering data from the nucleon [15]. In 
this work we analyze the impact of the off-shellness nature of the A in the EM responses for 
high momentum transfers by comparing the results obtained with different sets of parameters 
fitted to nucleon data. In this context we study, in particular, the differences with the 
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Figure 1: The Feynman diagrams corresponding to the two-body electroexcitation of the A. 

traditional Peccei Lagrangian approach, which corresponds to specific choices of the off- 
shell parameters and coupling constants. Recently attempts have been made to design new 
A interaction lagrangians "consistent" with the number of spin degrees of freedom of the 
A from a rigorously field theoretical point of view [13, 16, 17]. However analyses of pion 
photoproduction with these interactions have not been performed yet, and these studies 
should be done before attempting to implement them into the MEC operators. 

The structure of the work is the following: in sect. 2 we present the relativistic model for 
the A current and provide expressions for the ph matrix elements in the RFC In sect. 3 we 
present the results of the calculation of the nuclear response functions including the A current 
for several values of the momentum transfer, up to q = 3 GeV/c. We make contact with the 
scaling properties of the response functions and compare the results with those obtained in 
non- relativistic approaches for low and intermediate momentum transfer. Finally in sect. 4 
we present our conclusions. In appendix A we discuss the properties of the A propagator 
and provide a new method to derive its general form fulfilling the point transformation. In 
Appendix B we give details on the non-relativistic reduction of the A current. 



2 The relativistic A-isobar current 

In line with refs. [1, 2] we perform our analysis in first-order perturbation theory, namely we 
consider the contributions arising from diagrams with only one pionic line — our conventions 
are discussed at length in the references cited. We then evaluate the contribution of a virtual 
A to the longitudinal and transverse response functions 
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Figure 2: The particle- hole Feynman diagrams corresponding to the two-body electroexci- 
tation of the A. 



linked to the hadronic tensor W^ v '. In the symmetric (Z = N) RFG model the one-body-A 
interference contribution to the hadronic tensor reads 1 
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where (E p , p) = P M = (H + QY = (E h + u, h + q) and the lower limit of the integral, h , is 
the minimum momentum required for a nucleon to participate in the process (see refs. [2, 18] 
for details and explicit expressions). The ph one-body EM current is given by 



# B (p, h) = u(p) (Frf + i^-a^Q^j «(h) , 



(4) 



where F 1 and F 2 are the Dirac and Pauli form factors and u(h) = u(h, s^, ^) the free Dirac 
spinor. The current j A (p,h), associated with the diagrams of fig. 2, is linked to the ph 
matrix element of the A current operator, j A (Q) through 
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The sum J2k<k F becomes, in the thermodynamic limit, an integral over the momentum in 
the range < k < kp, and over the angular variables 9k,(f>k- The first and second terms in 
eq. (6) represent the direct and exchange contribution to the matrix element, respectively. 
We may write the 7./VA Lagrangian in the general form [15] 
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1 We assume q > 2kp, where &p is the Fermi momentum, so that no Pauli blocking is present. 
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4ata = j^^r Q a ,(z 2 , A)^(d v N)F^ + h.c. (9) 



C?,na = -j^^rQa^A^TiNd^ + h.c. , (10) 

where ip a is the A-field, iV the nucleon field and F Ufl the EM field tensor. We denote by TJ 
the | — > | isospin transition operators [19] for a = 1, 2, 3. The Lagrangians (8-10) coincide 
with the expressions given in refs. [5, 14, 20]. On the contrary, the term in eq. (9) differs 
from ref. [15] in a global sign. As it will be shown later, this difference makes a negligible 
effect on the global contribution to the MEC(A) transverse response, but it causes a very 
significant reduction in the contribution to the longitudinal response. 
The tensor Q^ u (z, A) can be written in the general form 

0^(z,A)=g^ + z + ^(l+Az)A llllv , (11) 

where z is the so-called off-shell parameter and A is an arbitrary parameter related to the 
"contact" invariance of the Lagrangian [5, 15]. 

To derive the A current we need also the 7riVA Lagrangian, given by 

CirNA = — ^e^(^, A)d»7T ■ T^N , (12) 

where it is the isovector pion field, and we have introduced an off-shell parameter z n for the 
7riVA vertex. Finally, as in ref. [2], we use pseudo- vector coupling for the nNN vertex. 

The corresponding two-body Delta current is obtained by computing the S'-matrix ele- 
ment for the elementary virtual photo-absorption process by two nucleons N± + N 2 + 7 — > 
N[ + N2 (see fig. 1). The corresponding current function can be written as 2 (see ref. [2] for 
a detailed definition of relativistic currents): 

f A (p'i, P' 2 , Pi, P2) = f -^G,{K 2 )K 2a Q v u{ v ' 1 ) [X^(P[, P x ) - X^(P[, P,)} u( Pl ) 

xu(p' 2 ) 75 F 2 T a u(p 2 ) + (1 <— > 2) , (13) 

where we use the Einstein convention for the Lorentz indices and for a sum over a repeated 
isospin index a = 1,2,3. Moreover, = P/ M — P- 1 (with i = 1,2) are the pionic four- 
momenta (see fig. 1) and 

GAK) = (14) 



K 2 — m 



is the propagator of a pion carrying four- momentum . The tensor X^ u is defined as 
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2 Here we use the Bjorken and Drell conventions [21], whereas different conventions were used in [22]. 
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The isospin sums in eq. (13) can be performed using the relations 
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In eq. (15) the three amplitudes contributing to the N-A vertex [7, 15] are taken into account 
and the values of the off-shell parameters z n , z±, z 2 and z 3 will be discussed later. 

In what follows we start by considering the A as a stable particle with mass toa, and we 
will later include a width to account for its decay probability in the resonance region. The 
isobar propagator can be expressed in general as a sum of two terms 
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is the usual Rarita-Schwinger (RS) propagator tensor and 
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the piece of the propagator that depends on the parameter A. Note that the global sign 
in the RS term differs from the expressions given in [5, 14, 15]. This can be due just to 
a different choice of phase in the definition of the propagator. As shown in Appendix A, 
the choice of phase in this work coincides with the one that provides the standard form of 
the nucleon propagator, Sn(P) = 1/(P — m N ). In appendix A eq. (20) is derived using a 
contact transformation of the RS propagator 3 . Moreover, because of the contact invariance, 
the physical properties of the field can be shown not to depend on A. As a further test of 
our calculation, we have checked that the results for the T response do not depend on A. 
Hence in what follows we fix A = —1 so that the complete A propagator is simply reduced to 
the Rarita-Schwinger expression and therefore omit the explicit A-dependence in the tensor 

It is immediate to check that the current in eq. (13) is conserved, 



Q*Ja(p1»P2>Pi»P2) = °> 



(21) 



so that the A current does not affect the gauge invariance of the RFG model. 

3 Note that there is an error in the relative sign between the and G^ p pieces given in refs. [5, 12]. 
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The particle-hole matrix-element diagrams are obtained from eq. (13) by setting 
{Pi, s[, t[) = (P, s p , t p ) , (P u Sl , tj) = (H, s h , t h ) 

(P2,s 2 ,h) = (P 2 ! ,s' 2 ,t 2 ) = (K,s k ,t k ) (22) 
for the direct term and 

(Pi, si, t[) = (P, s p , t p ) , (P 2 , s 2 , t 2 ) = (H, s h , t h ) 

(Pusuh) = (P 2 \,s' 2 ,t' 2 ) = (K,s k ,t k ) (23) 

for the exchange term, and by summing over the hole momentum k, spin s k and isospin t k . 
The isospin trace yields a vanishing direct matrix element, since Tr r a = TrT^T a = 0. In the 
exchange channel, recalling eq. (16), we get 

^2xt p T a T3XtkXt k T aXt h = ^x\ k T?,Tlxt h x\ v T aXt k = -Xt p r zXt h (24) 
tk t k 

and 

x\TzTixt k x\jaXt h = Xt k T a Tlxt h Xt p TaXt h = . (25) 
t k tk 

As a consequence only the two "vertex" diagrams (a) and (b) in fig. 2, contribute to the 
process, whereas the "self-energy" diagrams (c) and (d) give a vanishing contribution to the 
responses. 

Finally, by taking the thermodynamic limit, we get for the ph current function in eq. (5): 
/ttJVa/4 t r d 3 k m N ^ 
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correspond to the forward- and backward-going diagrams (a) and (b) of fig. 2, respectively. 

The decay of the A-resonance into a physical N-n state should be taken into account 
above threshold, i.e., P A > ( m ?r + m N) 2 - This is accounted for by modifying the propagator 
in eq. (19) to include a finite width according to the following prescription [6] 
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is the energy-dependent width. In the above p* is the momentum of the final pion resulting 
from the A decay (in the A-system) and p r ^ s is its value at resonance. Moreover we take 
r = 120 MeV. 



3 Results 

In this section we show the contribution of the A current to the longitudinal and transverse 
quasielastic response functions in the lp-lh channel. 

3.1 Response functions and off-shell dependence 

The following coupling constants have been used in the calculation [6, 15]: 

/ = V4tt x 0.08 , f nNA = 4 x 0.564 , G ± = 4.2 , G 2 = 4 , G 3 = 1 . (31) 

The value of G 3 is arbitrary owing to the lack of experimental information (see below). 
Moreover, although not explicitly indicated in the above formulae, the following monopole 
form factor 

A 2 _ m 2 

F nNN (K) = F nNA (K) = A2 _£ 2 , (32) 

with A = 1300 MeV, has been used. For the single-nucleon current we have for simplicity 
adopted the Galster form factor parameterization [23]. The Fermi momentum is chosen to 
be hp = 237 MeV/c, namely, representative of a typical sd-shell nucleus. 

In fig. 3 we plot the one-body-A interference contributions to the longitudinal and trans- 
verse responses as functions of the energy transfer u for various values of the momentum 
transfer q, ranging from 0.5 to 3 GeV/c. The separate contributions of the three Lagrangians 
in eqs. (8,9,10) are displayed. 

Regarding the longitudinal response (left panels), the G\ and G 2 pieces are similar in 
magnitude and tend to cancel for high q, whereas the G 3 term is negligible. Note that this 
result is very different from the one obtained with the Peccei Lagrangian (corresponding to 
the G\ term only). However, as will be shown later, the whole contribution to the response 
is very small. 
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Figure 3: The contribution of the A current to the longitudinal (left panels) and transverse 
(right panels) responses plotted versus 00. Here kp=237 MeV/c. The separate contributions 
of the first (Gi, solid) second (G 2 , dashed) and third (G 3 , dotted) terms of the current are 
displayed. The off-shell parameters are taken as Z\ — z 2 — z 3 — z^ — — 1/4. 
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Figure 4: Global longitudinal and transverse responses plotted versus u. Solid: RFG trans- 
verse response; dashed: RFG+MEC(A) transverse response; dotted: RFG longitudinal re- 
sponse; dot-dashed: RFG+MEC(A) longitudinal response. 

In the transverse channel (right panels) the G\ term clearly dominates, although at high 
q the contribution of G2 becomes significant and tends to cancel the first one. Hence in this 
case the results almost coincide, for q not too high, with the results obtained with the Peccei 
Lagrangian. 

To clarify the effects introduced by the A current, we show in fig. 4 the "total", namely 
RFG+MEC(A), longitudinal and transverse responses compared with the pure RFG ones. 
The same values of q as in the previous figure have been considered. We note that the con- 
tribution of the MEC(A) to the longitudinal response is always negligible. In the transverse 
channel the contribution of the A is larger and negative, that is, the interference between the 
one- and two-body currents matrix elements is destructive and therefore reduces the total 
from the purely nucleonic answer. At low momentum transfers the fractional contribution 
arising from MEC(A) contributions is relatively large and then it slowly decreases as the 
momentum transfer increases into the several GeV regime. Specifically, at the selected values 
of q = 0.5, 1, 2 and 3 GeV/c the net effect of the MEC(A) contributions is 19%, 18%, 10% 
and 4%, respectively. 

It is also interesting to compare the role of the MEC(A) with the other pionic MEC 
(seagull and pion-in-flight), calculated in the same relativistic model in [1, 2]. As shown 
in fig. 5 the A contribution is larger than the pionic one for lower q, where the total MEC 
are sizable, although for higher momentum transfers the difference in magnitude between 
A and pionic MEC contributions decreases and the two contributions tend to cancel. As a 
consequence, the impact of the total MEC at high q (say q > 2 GeV/c) is almost vanishing 
in the lp-lh sector. 
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Figure 5: MEC(A) contribution to the transverse response (solid) compared with the pionic 
seagull plus pion-in-flight contribution (dashed). 

We now investigate the effect of changing the off-shell parameters Z{. The values Zi = 
—0.25 used in the above results are unnecessarily restrictive, since they were fixed by Peccei 
under the condition 7^6^(2, A) = 0, which does not correspond to the most general form of 
Q M j, consistent with the point transformation. Different ranges or sets of the parameters z n , 
z\ and Z2 have been obtained by fitting pion photo-production data on the nucleon [14, 24] , 
while the value of Z3 needs electrons to be fixed. The uncertainty in these parameters arises 
from the treatment of the unitarity constraint, due to both theoretical ambiguities and to 
lack of precise experimental data. 

Since a unique determination of the off-shell parameters is not available, we have studied 
the dependence of our results upon a variation of them. In fig. 6 we present the A contribution 
to the longitudinal and transverse responses for q = 1 GeV/c using four different choices of 
the parameters (z n , z±, z%, 23): 

a) (-1/4,-1/4,-1/4,-1/4) , b) (-1/4,0.1,2.25,-1/4) , 
c) (-1/4,0.3,2.25,-1/4) , d) (-1/2,-1/2,-1/2,-1/2) . (33) 

Set a), used in the previous figures, reduces to the standard Peccei Lagrangian in the case 
G2 = G3 = 0. Sets b) and c) are determined in [15] by fitting Compton scattering cross 
sections on the nucleon (note that in this case the G3 term does not contribute) and set 
d) yields the usual non-relativistic 7^ A vertex, namely 6^ = g^. Although several other 
choices of parameters have been suggested in the literature [14], we believe that the results of 
fig. 6 give an indication of the uncertainty in the response functions associated with different 
off-shell prescriptions. It appears from fig. 6 that the fractional uncertainty is much larger 
in the longitudinal channel, where, however, the contribution to the response is negligible. 



10 



> 

05 



> 




-0.15 



250 300 350 400 450 500 550 600 




250 300 350 400 450 500 550 600 



U3 [MeV] 



Figure 6: MEC(A) contribution to the longitudinal (upper panel) and and transverse (lower 
panel) responses plotted versus u. The momentum transfer is q — 1 GeV/c and results are 
shown for the four sets of the off-shell parameters discussed in the text. 
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On the other hand the transverse response, for which the MEC(A) contributions are more 
important, is rather insensitive to the different choices (the effect being at most of the order 
of 6-7%). 

In [17] a correspondence between classes of "consistent" and "inconsistent" lagrangians 
was found by a redefinition of the spin-3/2 field. A new contact interaction that does not 
involve the spin-3/2 field appears. In our case this would mean that a new MEC term 
of contact type, dependent on the off-shell parameters, would arise and should be added 
to the "seagull" current, cancelling the dependence of the total responses on the off-shell 
parameters. 

In the next subsection we shortly address the scaling behavior of the MEC(A) responses. 



3.2 Scaling 

The scaling phenomenon has been presented in detail in refs. [25, 26, 27]. Here we just recall 
the basic definitions which are of use for the discussion that follows. We only focus on the 
transverse channel since, as previously shown, the MEC(A) are irrelevant in the longitudinal 
channel. 

Scaling of first kind occurs if the scaling function 

f {q ^ ) = kF Wi^) (34) 

becomes a function of one single variable, the scaling variable, and independent of q. Such 
behaviour is known to occur for large q in the region below the QEP. In eq. (34) G T (q, uj) is 
the relevant single-nucleon EM function (see ref. [27] for its explicit expression). 

Several different scaling variables exist in the literature, all of them coalescing into one 
- or being simply related to each other — for high enough momentum transfers. In the 
quasielastic peak region the natural scaling variable turns out to be [28, 29] 

(35) 

where T = \Jh$ + m 2 N — m N is the minimum kinetic energy required to a nucleon to take 
part in the process. The +(— ) sign in eq. (35) refers to the right (left) of the quasielastic 
peak. The analysis of the World data [25, 26, 27] shows that scaling of first kind is reasonably 
good for ij} < and badly violated for > 0. 

Scaling of second kind corresponds to the independence of the function f T on the specific 
nucleus, namely on the Fermi momentum. The analysis of the existing data points to an 
excellent fulfillment of this scaling in the region ip < and to a not very dramatic breaking 
of it for ip > 0. When the two kinds of scaling occur the response is said to "superscale" . 

The relativistic Fermi gas model fulfills both kinds of scaling, by construction, yielding the 
scaling function / = 3(1 — ip 2 )/4:. The observed superscaling behavior of the experimental 
data [25, 27] offers a clear constraint on the size allowed for nuclear correlations and MEC 
contributions, since these may break the scaling behaviour (of both kinds), and can therefore 
be used as a test of the reliability of the model. It is then natural to explore the scaling 
behaviour predicted in the present model. 
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Figure 7: The scaling function f T plotted versus ip at kF = 237 MeV/c for various values of 
q (upper panel) and at q — 1 GeV/c for various values of hp (lower panel). The function f 
(thick solid line) refers to the free Fermi gas scaling function. 

In [1] the evolution with q and hp of the pionic MEC and correlations has been already 
explored in detail: it has been proven that these contributions satisfy, at momentum transfers 
above 1 GeV/c, scaling of the first kind, but that they violate the second-kind scaling by 
roughly three powers of hp, which is strong enough to be seen in existing data. However, the 
size of the scale breaking predicted in [1] in the scaling region (below the QEP) was small 
enough to be compatible with the high quality high-g data. 

In fig. 7 we display the scaling function f T versus ip for different values of q (upper panel) 
and hp (lower panel). It appears that both kinds of scaling are violated by the MEC(A) 
contribution, something not evident in the high quality World data, at least within exist- 
ing experimental uncertainties. When the "total" lp-lh response is formed a quantitative 
analysis shows that the MEC(A) contributions in the transverse response play some role, 
but that overall first-kind scaling is quite good for — 1 < ip < —0.5 and q < 1 — 2 GeV/c 
and is violated at a 5-10% level at the QEP. On the other hand, the second-kind scaling 
violations from this lp-lh model can be significant for all values of vp, since they modify the 
scaling function (that scales, by construction) by three powers of kp. The figure shows a 
generous range of Fermi momenta; these correspond to going from 4 He at kp = 200 MeV/c 
to very heavy nuclei at kp = 250 MeV/c and then well beyond to 300 MeV/c to explore 
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more fully the second-kind scaling behaviour. If we restrict our attention to nuclear species 
where high quality data exist then we should compare the kp = 200 MeV/c results with 
those for kp = 250 MeV/c and the lower panel in the figure shows that violations of about 
14% are predicted at the QEP (actually using modeling of this type for helium is probably 
stretching things somewhat: if the comparison is made for "real nuclei" such as 12 C ver- 
sus 197 Au for which excellent data exist then an effect of roughly 8% is predicted). While 
second-kind scale breaking of this magnitude should be visible in the data one should be 
careful before drawing premature conclusions, since one knows that additional scale-breaking 
contributions arise from mechanisms outside the present model. In particular we know from 
past work [6, 22, 30] that 2p-2h MEC contributions also typically lead to effects that go as 
kp, but these add and therefore tend to cancel the contributions discussed in the present 
work. The net effect, in the scaling region at least, is expected to amount to a few percent 
and therefore be consistent with the relevant World data. 

It appears that the scaling behaviour of the MEC(A) contributions is different from that 
of the single-nucleon RFG response — namely, if the scaling functions and variables are 
defined as in past work to make the latter scale, then the two-body MEC contributions in 
general do not. First-kind scaling violations (see the upper panel in fig. 7) arise because the 
one- and two-body currents do not in general have the same (q, u) dependences, and only 
one choice can be made when defining the scaling function. Namely, in standard treatments 
it is the one-body (single-nucleon) cross section that is divided out to produce scaling for 
the dominant impulse approximation contributions. Similarly, second-kind scaling violations 
arise when the density dependences of the various contributions are different. 

The high-quality World data present a mixed picture: for negative values of ip both 
kinds of scaling are reasonably respected, although scaling of the second kind appears to 
be better than scaling of the first kind in this region. For ip positive, scaling of the first 
kind becomes quite bad, partially because pion production including via the delta and other 
baryon resonances becomes important and, for the reasons stated above, has a different 
(q,u) dependence. Additionally, the MEC(A) effects under study in the present work can 
contribute to this first-kind scale breaking. Were the latter to be the only contributions to 
be present along with the one-body responses, then we would conclude from this study that 
overall first-kind scaling is quite good for — 1 < ip < —0.5 and q < 1-2 GeV/c, but could 
be violated by as much as 5-10% at the QEP. However, the full assessment of which contri- 
butions produce the observed result cannot be made before all effects are included, namely 
from impulse approximation nucleonic currents, lp-lh MEC effects from pion-exchange con- 
tributions including those involving the A (this work), 2p-2h MEC effects (work in progress) 
and meson production. 

Similar statements can be made with respect to second-kind scaling violations. The 
MEC (A) contributions studied in the present work (see lower panel in fig. 7) can be signifi- 
cant for all values of ip, since they modify the scaling function (that scales, by construction) 
by three powers of kp. 
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Figure 8: MEC(A) transverse response computed for several values of q and hp. Solid: exact 
relativistic result, dashed: non-relativistic result. Static propagators without a 7r iViV form 
factor have been used. 
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3.3 Non-relativistic reduction 



To finish with the discussion of the results, in this section we analyze the effects associated 
with different types of non-relativistic approaches. Focusing on the MEC(A) contribution 
to the hadronic (e, e') observables one deals with the single-nucleon EM current operator 
and the two-body A current. Improved "semi-relativistic" (SR) versions of the former have 
been presented in [31, 32, 33] by expanding only in the dimensionless parameter r\ = p/m,N, 
being p the struck nucleon three-momentum. By this method a SR form of the transverse 
OB current can be written in terms of the non-relativistic one as 

&h) = ifflp,h), (36) 
Vl + r 

where r = \Q 2 \/4m 2 N . A similar expansion has been carried out for the pion-in-flight and 
contact MEC in ref. [34] and for the one-body A electroproduction current in [35]. In the 
case of the MEC, a semi-relativistic expression similar to eq. (36) was proposed and tested 
in [2] 

j^ C (p,h) = _^j^(p,h). (37) 

Although these SR MEC currents compare better than the NR ones to the exact relativistic 
calculation, they are not as good as that the SR OB current since, as illustrated in [2], an 
additional /c^-dependent normalization factor N(q,u, hp), arising from the integration over 
the Fermi sea, should be present in eq. (37). Nevertheless the simplicity of these SR currents 
make them suitable for easy incorporation in existing non-relativistic models. 

The relativizing factor (1 + t) -1 / 2 takes care of spinology and normalization properties 
that cannot be neglected in any relativistic model. Therefore, consistently with eq. (37), we 
propose a similar expression for the semi-relativistic A current 

jf*(P> h ) = - J=^*(P> h ) - (38) 

Vl + T 

where j^ R (p, h) is the standard non-relativistic A current commonly considered in the liter- 
ature [36, 37], usually obtained from NR sets of 'jNA and ttNA Lagrangians. In appendix 
B we derive its expression by performing a non-relativistic reduction of the relativistic A 
current considered in this work. 

To illustrate clearly the impact of relativity we show in fig. 8 the MEC(A) transverse 
response corresponding to the fully relativistic calculation (solid lines) and to the standard 
NR reduction (dashed lines) using the Fermi gas model of ref. [38]. In both cases no irNN 
form factor has been used and the static limits of the pion and A propagators have been 
assumed. We observe that the two calculations give the same results for small density and 
momentum transfer and start to differ as q and hp increase. In particular the results of fig. 8 
provide a test of the multi-dimensional numerical integration procedure used in this work, 
since the A integrals over k in the non-relativistic calculation of [38] are analytical. Moreover 
the agreement of both results for low momenta also represents a test of the applicability of 
the non-relativistic A current in this kinematic domain. This result is not obvious a priori 



16 




Figure 9: MEC(A) transverse response for q — 1 GeV/c. The result corresponding to 
the fully relativistic Fermi gas (RFG) is compared with the non-relativistic (NR), including 
relativistic kinematics (RK), semi- relativistic approach for the one-body current (SROB) 
and for the A-current as well (SR) (see text for details). 

due to the somewhat coarse approximations assumed in the derivation of the NR current 
(see Appendix B). 

The shrinking of the response domain in fig. 8 with increasing q arises from the relativistic 
kinematics in the energy-conserving delta function appearing in the responses. As in past 
work this effect can be accounted for approximately by the replacement A — > A(l + A), with 
A = u/2rriN, in the non-relativistic calculation. 

This is shown in fig. 9, where we select the case q = 500 MeV/c and compare the NR 
calculation (dashed line) with the same calculation but using relativistic kinematics (RK) 
(dotted line). Apart from this effect, the enhancement of the non-relativistic calculations 
with respect to the relativistic one is a consequence of the reduction of the single-nucleon and 
A current operators. The semi-relativistic approach for the single-nucleon current (SROB) 
(dot-dashed) improves the agreement with the fully relativistic result. Similarly the response 
obtained by using the SR current (38) (thin solid line) gets a little bit closer to the relativistic 
result, although a difference between them still exists. Better agreement between the exact 
and the relativized models for the MEC(A) responses are obtained in the limit hp — > in 
the quasielastic peak region. 

4 Conclusions 

In this paper we have studied the role of the MEC(A) two-body current in lp-lh quasielastic 
electron scattering. A fully relativistic, gauge invariant, analysis has been performed in the 
context of the RFG. The most general Lagrangian has been considered, allowing for different 



17 



off-shell parameters in each one of the vertices, and the roles played by its various pieces have 
been analyzed. The results have been compared with various non-relativistic reductions and 
the behavior of the A contribution with respect to the scaling phenomenon has been also 
explored. 

In summary, our results show that: 

1. The A contribution is only important in the transverse channel, where it represents 
the dominant correction to the free Fermi gas, at least for q < 1-2 GeV/c. 

2. In the transverse channel the main contribution comes from the first term of the La- 
grangian (Peccei). This is strictly valid for not too high values of q. 

3. In the longitudinal channel the global A contribution is negligible. In this case, how- 
ever, the first and second terms of the Lagrangian tend to cancel. 

4. The transverse response is quite insensitive to the different sets of off-shell parameters 
used, whereas the longitudinal one shows a much more pronounced sensitivity to them. 

5. The MEC(A) contributions by themselves clearly break both kinds of scaling be- 
haviour. When combined with the RFG the total scaling breaking is on the edge 
of being an effect that should be seen in the high quality high-g data; however, it 
is known that contributions from 2p-2h MEC tend to cancel the lp-lh MEC effects 
studied here, and thus the net scale breaking is in reasonable accord with the data. 

6. The effects of relativity are already sizable at q — 500 MeV/c and they can only 
be accounted for approximately by using relativistic kinematics and correcting the 
operators with a semi-relativistic prescription similar to the one found in the OB 
current and MEC operators. 

Finally, with the present relativistic model we have completed the gauge-invariant MEC 
model of ref. [2] to a fully consistent, gauge invariant approach to the electro-nuclear re- 
sponses in the one-nucleon emission channel, including all the diagrams with one-pion ex- 
change and intermediate excitations of A isobar, which can be applied for high-|Q 2 | values 
in the region of the quasielastic peak. 

A A-propagator and contact invariance 

In this Appendix we derive the general relativistic expression for the A propagator, discussing 
in particular its sign and the invariance of the theory under contact transformations. 

The current amplitudes of fig. 1 are computed using the standard Feynman rules which 
are essentially based on contraction of pairs of fields or propagators. The A-propagator 
in coordinate space is then derived, in analogy with the Feynman propagator for spin- 1/2 
fermions. It reads 




(39) 
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where 



= E 



P,SA 



c SA u a (p, s A )e~ w - x + 4 A v a (p,s A )e 



+iP-X 



(40) 



is the free massive spin-3/2 Rarita-Schwinger field, resulting from the coupling of a spin- 1/2 
and a spin-1 object: 



u a (p, s A ) = E(| slA H s A)e a (p, A)«(p, s) , 

X,s 



(41) 



■u(p, s) being a Dirac spinor of mass m A and e a (p, A) the basis vectors. 

The Rarita-Schwinger propagator in momentum space is the 4-dimensional Fourier trans- 
form of G R s(X) 



G 

and it reads 



-iPX 



(2nY 

P + m A 
P 2 — m A 



■ _ i _ 2 P*P?_ _ 1 - 
3 3 m A 3 mA 



(42) 



(43) 



Note that the minus sign in eq. (43), directly arising from the definition in eq. (39), is 
opposite to the one in refs. [5, 14, 15], whereas it agrees with ref. [39]. A simple check of the 
sign can be done by evaluating the right-hand side of eq. (39) for t > t'. Inserting the field 
expansion of eq. (40) in eq. (39) we obtain 



(X-X') = J ^^E^(P^A)^( P)S A)e-^^'). 



(44) 



The spin sum inside the integral reads 
E m "(P' s a)m /3 (p, s a ) 

EEE(| s ^l| s A)(^ / lA'|| S A)e a (p,A)e' 3 (p,A')MP ! «P^') 



(45) 



sa As X's' 



This sum is particularly simple in the static limit p = and for the components a = (3 = 3: 
in this case e a (p, A) and e /3 (p, A') vanish unless A = A' = 0. Moreover the Clebsch-Gordan 
coefficients select s = s' = s A , yielding 



5> 3 (0,s a )H 3 (0,sa) = E(i sl0 H s > 2 u(0,s)u(p,s) 



SA 



- Y u(0, s)u(0, s) = , 

3Y 3 2m A 



(46) 



where we have used 



L1 io|§±> = 



» 2 2 I 2 2 



, and where the positive-energy 



projector for spin 1/2 and mass m A is meant to be evaluated at p = 0. The same result is 
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obtained from the RS propagator in eq. (43) setting p = and a — f3 — 3, performing the 
contour integral over p in eq. (42): 

f dp f> + m A imt f + m A 

J 27rP 2 -m 2 A + te " 2E P { ] 



and using 

2 f + mi 



r 33 rp^ - ^ + mA 



y 33 -^(7 3 ) 2 



3P 2 



(48) 



Let us now discuss the so-called "contact-invariance" , induced by the operator 



0»„(b) = g\ + b -S\ , (49) 

where b is an arbitrary parameter and S^ u = 7^7^ is the generator of the transformation, 
fulfilling the identity S 2 = AS. 

A contact transformation of the RS field is defined by 

vfj' a = O a p {b)^ . (50) 

Some useful properties of the contact operators are: 

i) 0(a)0(b) = 0(a + b + 2ab) (51) 

a > "< 6 » = °(-rT2s) < 52 > 

iii) 0(a)0- 1 (b)=o(j^\ . (53) 

The contact invariance follows from the fact that the A-dependent free Lagrangian of the 
|-fleld becomes independent of the parameter A, after a redefinition of the RS field [15]: 

A" = O a p {A)^ . (54) 

Therefore, in terms of the contact transformation eq. (50) we can write the new field in 
matrix form as: 

A = 0(A)0-\b)il)' (55) 

and using eq. (53) we obtain that contact invariance is fulfilled if the parameter A changes 
to A' = (A — b)/(l + 2b). 

We now derive the general form of the A-propagator G(A) for an arbitrary parameter A 
in terms of the usual RS propagator in eq. (43). From the definition in eq. (39) (for simplicity 
we omit the X, X' variables) we get 

iG aP {A) = (0\T{o£(A)ATO;j(A)}\0) 
= O£(A){0\T{A*Sr}\0)O;}(A) 

= O-l(A)iG^(0)O^(A), (56) 
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where G^(0) is the propagator for A = 0, obtained in terms of the free Lagrangian for the 
A a field. The above equation can be written in matrix form introducing the 4x4 matrix of 
the propagator G(A) = G^ U (A) 

G(A) = O-\A)G(0)O-\A) . (57) 

The RS propagator in eq. (43) corresponds to A = — 1, which also can be related with the 
A = case 

G(-l) = O-\-l)G(0)O-\-l) . (58) 
Using the two above equations we can write 

G(A) = 0- 1 (A)0(-l)G(-l)0(-l)0- 1 (A) 

= 0(B)G(-1)0(B) , (59) 

where 

B = (60) 
1 + 2A y J 

and use has been made of eq. (53). 

From eq. (59) it is immediate to show that the combination X(A) = Q(z, A)G(A)Q(z' , A), 

appearing in the MEC(A) current, is independent of A. We note that the operator in eq. (11) 

can be expressed, by means of eq. (51), as the product of two contact operators 

Q{z, A) = 0{2z)0{A) = 0(A)0(2z) . (61) 

Hence 

X(A) = 0(2z)0(A)G(A)0(A)0(2z') = O(2z)G(0)O(2z) = X(0) . (62) 

To obtain the explicit form of G(A), we note that the Rarita-Schwinger propagator G(—l), 
given by eq. (43), can be more conveniently written as 



G(-l)= f + mA 



1--S — \u —(V - W) 

3 3mi 3mA 



(63) 



P 2 — m\ 
having introduced the matrices 

U» v = P^P U , V\ = i^P u , W» v = P^i u . (64) 
From eq. (59) it follows that 

G(A) = G(-l) + | [SG(-l) + G(-1)S] + ^SG(-1)S = G RS + G A . (65) 
After straightforward 7-matrix algebra the following useful relations can be deduced: 

SG(-l) + G(-1)S = -^(V + W -m A S) (66) 

SG(-1)S = -^ T (2W- PS -2m A S) (67) 

OTfl A 

which finally yield 



Ga = - — 2~B 
3m A 



V + W- m A S + l -B(2W- PS - 2m A S) 



(68) 



This, exploiting eq. (60), coincides with eq. (20). 
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B Non-relativistic reduction of the A-current 

In this appendix we derive in detail the non-relativistic reduction of the A-exchange current 
corresponding to the Peccei-like 7./VA vertex. 

Let us consider the ^NA Lagrangian as given by eq. (8). The corresponding two-body 
current can be written as follows 

#(p'i, p' 2 , Pi, P 2 ) = ^ f -^-GAK 2 ) [A»T a Tt + STsTj] u(p 2 ) 7 5 fiVutfr) 

+ (1 <=► 2) , (69) 
where we have introduced 

A" = ^p[)K 2a e^G^(P 1 + Q)Q u (& p ^ u -& pu iniMPi) (70) 
B» = u(p[)K 2al5 QAl u e pp -r&nG^(Pi-Q)&^u(Pi)- (71) 

To compare with the Peccei vertex the tensor G Q/3 in eq. (11) should be taken for z = —1/4 
and the A propagator is the Rarita-Schwinger expression in eq. (19). 

In what follows we invoke the static limit as usually considered in standard non-relativistic 
calculations. In this case the spin dependence in the tensor a/3 is neglected, i.e., <d af3 = g af3 . 
Moreover, the pion propagator simply reduces to 



G - {K) = (72) 

while the A propagator only contributes for space indices becoming 

Gtj(P) = Uj + , ij = 1, 2, 3 . (73) 

J rriN — rriA V o / 

Assuming ~ (0, k 2 ) and ~ (0, q) (valid in the static limit), the space components of 
the four-vectors A^ and in eqs. (70,71) can be written in terms of space components only 



A 1 * «(p / 1 )^(5h + tt7^)^(-5V + ^V)75«(pi) (74) 

B l * Bfo'JfcfoMV + ^V) 76 U kl + \ lllk ) u(pi) • (75) 

m N — m A V 6 J 

Taking now the positive energy components of the gamma matrices, 

Ikli -> -cr k cri, 7i7 5 -> a h S k i + ^IHk -> t^ui ~ ^kim^m (76) 
we can write A\ B % in terms of matrix elements between the Pauli spinors, Xi> as follows 

A 1 = xi^Xi (77) 

m N — m A 

B* = xi^Xi, (78) 

ttin — rriA 
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where A 1 and B l are the non-relativistic reduced operators given by 
,2 2 i i 

2 2 2 i 

B l = —-k 2 iqjcrj + -k 2 jqj<T i + -e ikm k 2 kqjO'j(T m — ^jkm^kQj^i^m ■ (80) 

After some algebra involving vector relations and properties of the Pauli matrices, the re- 
duced non-relativistic A and B terms in their vector form result 
1 2 

A = -q x (k 2 x cr) - -iq x k 2 (81) 
— 1 2 

B = --q x (k 2 x cr) - -iq x k 2 . (82) 
Moreover, the matrix element of the second nucleon in the static limit reduces to 

w(p' 2 )75 ^2«(p 2 ) — ► u(p' 2 )k 2 ■ I1MV2) — ► X2 ( k 2 -o)x2- (83) 
Using the above results, the ith. component of the current can be written as 

^A = x'i t X 2 t ^AXiX2, (84) 
where, using eq. (16), 

-j7 -2 G\ f wN A f k 2 • cr (2) 1 r ( 2 ) r (1) (2) ] (1) 1 

J A ^ ^ 2,12 l 4r 3 k 2 ~ T X T °" X M X q 

+ (1—2). 

(85) 

This form coincides with the usual non-relativistic A current used in the literature [37]. 
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